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^ i Abstract 

We consider an extension of the standard model (SM) with charged singlet scalars and right 



Ph. 



handed (RH) neutrinos all at the electroweak scale. In this model, the neutrino masses are gener- 



Qh ated at three loops, which provide an explanation for their smallness, and the lightest RH neutrino, 

iVi, is a dark matter candidate. We find that for three generations of RH neutrinos, the model can 

be consistent with the neutrino oscillation data, lepton flavor violating processes, N\ can have a 

> 

I/"") ■ relic density in agreement with the recent Planck data, and the electroweak phase transition can be 

strongly first order. We also show that the charged scalars may enhance the branching ratio h — > 77 
■^j- ' , where as h — > ^Z get can get few percent suppression. We also discuss the phenomenological 

o' 

CO ■ implications of the RH neutrinos at the collider. 

PACS numbers: 95.35.+d; 98.80.-k; 12.15.-y; 11.30.Qc. 

Keywords: Charged scalars, WIMP's, Higgs diphoton decay, electroweak phase transition. 



Electronic address: ahrichc@univ-jijcl.dz 



^Electronic address: snasri@uae u.ac.ae1 



I. INTRODUCTION 

There are three concrete evidences for Physics beyond the standard model (SM): (i) non 
zero neutrino masses, (ii) the existence of dark matter (DM), and (iii) the observation of 
matter anti matter asymmetry of the universe. However, most of the SM extensions make 
no attempt to address these three puzzles within the same framework. For instance, in 
the minimal supersymmetric standard model (MSSM), the lightest supersymmetric particle 
(LSP) is a candidate for DM and, in principle, has the necessary ingredients to generate 
the baryon asymmetry of the universe (BAU), but it does not provide an explanation for 
why neutrino masses are tiny. Moreover, direct searches for supersymmetric particles have 
yielded null results so far. Another popular extension of the SM, is introducing very heavy 
right -handed (RH) neutrinos (mjv > 10 8 GeV, where small neutrino masses are generated 
via the see-saw mechanism |l| , and the BAU is produced via leptogenesis [2| . Unfortunately, 
such heavy particles decouple from the effective low energy theory and can not be tested at 
collider experiments. In addition, for m^ heavier than 10 7 GeV, the Dirac neutrino mass 
term induces large corrections to the Higgs mass, which can destabilize the electroweak 
vacuum 3|. 

Another possible way to understand the smallness of neutrino masses is to generate them 



l, 



radiatively. The famous example is the so-called Zee model |4j, where one augments the 
scalar sector of the SM with a Higgs doublet, and a charged field which transforms as a 
singlet under SU(2)l, which leads to non zero neutrino mass at one loop level. However, 
the solar mixing angle comes out to be close to maximal, which is excluded by the solar 
neutrino oscillation data [5|. This problem is circonvented in models where neutrinos are 
induced at two loops [6| or three loops |7H9j . One of the advantages of this class of models is 
that all the mass scales are in the TeV or sub-TeV range, which makes it possible for them 
to be tested at future colliders. 

In Ref. |7|], the SM was extended with two electrically charged SU(2)l singlet scalars and 
one RH neutrino field, N, where a Z2 symmetry was imposed to forbid the Dirac neutrino 
mass terms at tree level [7|. Once the electroweak symmetry is broken, neutrino masses are 
generated at three loops, naturally explaining why their masses are so tiny compared to the 
charged leptons as due to the high loop suppression. A consequence of the Z% symmetry 
and the field content of the model, N is Z 2 -odd, and thus guaranteed to be stable, which 



make it a good DM candidate. In Ref. [101 ]. the authors considered extending the fermion 
sector of the SM with two RH neutrinos, in order to for it to be consistent with the neutrino 
oscillation data, and they studied also its phenomenological implications. 

Here, we calculate the three loop neutrino masses exactly, as compared to the approximate 
expression derived in [7[ . We show that in order to satisfy the recent experimental bound on 
the lepton flavor violat ing (LFV) process such as /i — > ey ll| ; and the anomalous magnetic 
moment of the muon [12], one must have three generations of RH neutrinos. Taking into 
account the neutrino oscillation data and the LFV constraints, we show that the lightest 
RH neutrino can account for the DM abundance with masses lighter than 225 GeV. The 
presence of the charged scalars in this model will affect the Higgs decay process h — > 77 and 
can lead to an enhancement with respect to the SM, where as h — > 7Z is slightly reduced. 
In this model, we find that a strongly electroweak phase transition can be achieved with a 
Higgs mass of ~ 125 GeV as measured at the LHC [13l . Il4j . 

This paper is organized as follows. In the next section we present the model, and discuss 
the constraints from the LFV processes. In section III, we study the relic density of the 
lightest RH neutrino, and discuss the coannihilation effect due to the next lightest RH 
neutrino. Section IV is devoted to the study of the electroweak phase transition. In section 
V, we discuss the phenomenological implications of the RH neutrinos at electron-positron 
colliders. Finally we conclude in section VI. The formula of the three loop factor that enters 
in the expression of the neutrino masses is derived in Appendix A. In Appendix B, we give 
the shift in masses the gauge bosons and the scalars at finite temperature. 

II. NEUTRINO DATA AND FLAVOR VIOLATION CONSTRAINTS 

In this section, we define the filed content of the model, give the exact expression of the 
neutrino masses, and discuss the constraints from LFV processes. 

A. The Model 

Here we consider extending the SM with three right-handed neutrinos, TVj, and two elec- 
trically charged scalars, Sf and S 2 , that are singlet under SU{2)l gauge group. In addition, 
we impose a discrete Z 2 symmetry on the model, under which {S 2 ,Ni} — >■ {— S 2 , —N t }, and 



all other fields are even. The Lagrangian reads 

C = C S M + {f a pLlCiT 2 LpS+ + giaNiS^iaR 

+ iM^Af CW, + h.c} - V(®, S u S 2 ), (2.1) 

where, L a is the left-handed lepton doublet, C is the charge conjugation matrix, f a p 
are Yukawa couplings which are antisymmetric in the generation indices a and /3, and 
V(&, Si, S 2 ) is the tree-level scalar potential which is given by 

V($, S h2 ) = A (|$| 2 ) 2 - fi 2 |$| 2 + mlSlSi + m 2 2 S* 2 S 2 + XiS^Si |$| 2 + X 2 S* 2 S 2 1$| 2 

+ f (siSif + 1 (s*s 2 f + V i 2 s*SiS*s 2 + {\ s SiSiS;s; + h.c} . (2.2) 

Here $ denotes the SM Higgs doublet. It is worth mentioning that, the charge breaking 
minima are not possible due to the positive-definite values of A s and 7712; in addition to the 
conditions on the charged scalar masses m? Si = m 2 2 + Ai j 2'^ 2 /2 > 0. 

There are two immediate implications of the Z 2 symmetry imposed on the Lagrangian: 



First, if Ni is the lightest particle among N 2 ,N 3 , Si and 5*2, then it would be stable, 
and hence it would be a candidate for dark matter. Moreover, Ni will be pair produced 
and subsequently decay into N\ (or to N 2 and then to N±) and a pair (or two pairs) 
of charged leptons. We will discuss its phenomenology in section VI. 

The second implication, is that the Dirac neutrino mass term is forbidden at all levels 
of the perturbation theory, and Majorana neutrinos masses are generated radiatively 
at three-loops, as shown in Fig. [TJ 
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FIG. 1: The three-loop diagram that generates the neutrino mass. 



B. Neutrino mass 

The neutrino mass matrix elements arising from the three-loop diagram in Fig. [TJ are 
given by 

{M v ) aP = (t^Cl* f«iff>k9ij9kjF (m 2 Nj /m 2 S2 ,m 2 Si /m 2 S2 j , (2.3) 

where p, k(= e, p, r) are the charged leptons flavor indices, % = 1, 2, 3 denotes the three right- 
handed neutrinos, and the function F is a loop integral given in ( 1A8|) . which was approxi- 
mated to one in the original work J7|. Note that, unlike the conventional seesaw mechanism, 
the radiatively generated neutrino masses are directly proportional to the charged leptons 
and RH neutrino masses as shown in (I2.3P and (1A8|) . 

In general, the elements of the neutrino mass matrix can be written as 

(M v ) a = [U ■ diag(m 1 ,m 2 , m 3 ) ■ U T ] a(3 , (2.4) 



where U is the Pontecorvo-Maki-Nakawaga-Sakata (PMNS) mixing matrix 15[, which is 
parameterized in general by 
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(2.5) 
with Sij = sin(6ij) and Qj = cos(%), 5d is the Dirac phase; and a and /? are the Majorana 
phases. Using the experimental allowed values for s\ 2 = 0.320^o!oi7 5 s 2 3 = 0.43^o!o35 $i 3 = 
0.025±{j:oo3, lAm^l = 2.55±g;g| x 10~ 3 eV 2 and Am 2 21 = 7.62±° f 9 x 10" 5 eV 2 (lj, we can find 
the parameter space of the model that is consistent with the neutrino oscillation data. 

C. Experimental constraints 

Besides neutrino masses and mixing, the Lagrangian (12.1 p induces flavor violating pro- 
cesses such as £ a —¥ •yip if mg a > m^ , generated at one loop via the exchange of both extra 



charged scalars S 1 2 . The branching ratio of such process can be computed following [17] as 
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FIG. 2: The muon anomalous magnetic moment versus the (3f3ou decay effective Majorana (M u 
The blue lines represent their experimental upper bounds. 
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with k ^ ck,/3, a em is the fine structure constant and F%{x) = (1 — 6x + 3a; 2 + 2x 3 — 
6x 2 lna;)/6(l — x) 4 . For the case of £ a = lp = /i, this leads to a new contribution to the 
muon anomalous magnetic moment 5a ^ that is given by 
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(2.7) 
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In Fig. [21 we show a scattered plot of the muon anomalous magnetic moment versus 
the (3(3ou decay effective Majorana mass (M u ) ee . In our scan of the parameter space of 
the model, we took ms 1 2 > 100 GeV; and demanded that (12. 3p to be consistent with the 
neutrino oscillation data. From Fig. |2} one can see that most of the values of (M u ) ee that 
are consistent with the bound on 5a^ are lying in the range 10 -3 eV to ~eV. The current 
bound on (M u ) ee is approximately 0.35 eV 19| and it is expected that within few years a 



number of next generation f3f3o v experiments will be sensitive to (M u ) ee ~ 10 2 eV 20] . 



1 One has to mention that this results is different from Eq. (38) in [l8(, where the authors took the 
summation over the square of the gi a terms instead of the square of the their summation. The latter 
allows the parameter space of couplings to be enlarged. 
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FIG. 3: Different parameters combinations (as absolute values) that are relevant to the LFV con- 
strain on B(fi — y e^f), are shown where i2.3\) and \2.1$ are matched. 

Fig. E] gives an idea about the magnitude of the couplings that satisfy the constraints 
from LFV processes and the muon anomalous magnetic moment, and which also are con- 
sistent with the neutrino oscillation data. It is worth noting that when considering just 
two generations of RH neutrinos (i.e, gj, a = 0), we find that the bound on the experimental 
bound B (/j, — y ej) < 5.7 x 1CT 13 is violated 111] 2 . Therefore, having three RH neutrinos is 
necessary for it to be in agreement with the data from the bounds from LFV processes. 



III. DARK MATTER, COANNIHILATION EFFECT & INDIRECT DETECTION 



A. Relic density 



As we noted earlier, the lightest RH neutrino JVi is stable, and could be the DM candidate. 
In the case of hierarchical RH neutrino mass spectrum, we can safely neglect the effect of 
N2 and N% on N\ density. The Ni number density get depleted through the annihilation 
process N\Ni — y l a £p via the t-channel exchange of S^. For two incoming dark matter 
particles with momenta p\ and p2, and final states charged leptons with momenta k\ and 
k>2, the amplitude for this process is 
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(3.1) 



2 Although, we have considered the bound on B(t —y uj) < 4.5 x 10 8 12|, in our numerical scan, it docs 
constrain severely the parameter space of the model. 



where t = (pi — k\) 2 and u = (pi — k 2 ) 2 are the Mandelstam variables corresponding the t 
and u channels, respectively. After squaring, summing and averaging over the spin states, 
we find that in the non-relativistic limit, the total annihilation cross section is given by 

a NlNl v r ~ 2^ \9ia9ip\ — rT4 v n (3.2) 

a,p 487r(m| a +m^J 

with tv is the relative velocity between the annihilation iVi's. As the temperature of the 
universe drops below the freeze-out temperature Tf ~ m^/25, the annihilation rate becomes 
smaller than the expansion rate (the Hubble parameter) of the universe, and the iVi's start 
to decouple from the thermal bath. The relic density after the decoupling can be obtained 
by solving the Boltzmann equation, and it is approximately given by 

„ ,0 2x f x 1.1 x 10 9 GeV _1 , . 

s/g*M p i {cr NlNl v r ) 

1.28 XlQ- 2 . mNi x2 (l+m| /m^) 4 
~ V 1 \n n* |2 V135 GcvJ 1+mi /m 4 „ ' 

where < -u r 2 > ~ 6/x/ — 6/25 is the thermal average of the relative velocity squared of a pair 
of two iVi particles, M p i is planck mass; and g*{Tf) is the total number of effective massless 
degrees of freedom at Tf . 

In Fig. HI we plot the allowed mass range {m^^ms^} plane that give the observed dark 



matter relic density 21(. As seen in the figure, the neutrino experimental data combined 
with the relic density seems to prefer ms 1 > ms 2 for large space of parameters. However, 
the masses of both the DM and the charged scalar S 2 can not exceed m^ 1 < 225 GeV and 
ms 2 < 245 GeV, respectively. 

B. Coannihilation effect 

In computing the relic density in (13. 3p . we have assumed that there is a hierarchy between 
the three right-handed neutrino masses. However, if we consider the possibility for N 2 and/or 
N 3 being close in mass to N\, i.e A$ = {rn^ i —m^/m^ << 1, then coannihilation processes 
like N 1 N 2t 3 — > £ a £/3 might have important effects on the evolution of the iVVnumber density. 
The process NiSf- —> ^7 is suppressed by the large mass difference between S 2 and Ni 
and the smallness of the electromagnetic coupling compared to 0(g 2 ), and therefore its 
contribution to the coannihilation is negligible. 
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FIG. 4: The charged scalar masses ms 1 (red) and ms 2 (green) versus the lightest RH neutrino 
mass, where the consistency with the neutrino data, LFV constraints and the DM relic density 
have been imposed. 



Following 22( , the coannihilation effect could be included by re-evaluating the relic den- 
sity 03. 3p using the effective annihilation cross section and multiplicity 



°eff(x) = E- k f~ (1 + A,) 3/2 (1 + A*) 3/2 e-^+^a (N.N, -> i~£+) 



(3.4) 



^//(^) = E-2(l + A t ) 3 / 2 e-^, 

For freeze-out temperature much smaller than m^ 1 , the effective cross section can be written 
as a e ff(x)v r = a e ff(x) + b e ff(x)v 2 + O (v*), where 
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(3.5) 



Here the factors a^ and b^ correspond to the two first terms in the velocity expansion of 
a (NiNk — > i-a^p) v r (he, the s and p wave terms, respectively), given by 
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Thus, the coannihilation effect on the relic density could be accounted for by just multiplying 
the couplings term J2 a a \9iag*a\ 2 i n A3.3J) by the factor 

Here a/^ = m^/Ti is the freeze-out temperature defined using a e ff(x)v r instead of cr^^Vr, 
and the integral functions I a (x) and /& (x) are given by 

OO /"OO 



I a (x)—x a e ff(z)z dz, h(x) = 2x / b e ff(z)z dz, (3.8) 

J X J X 

which, in general, causes a shift in the freeze-out temperature. In our case, we find that Xf 
gets lowered by less then 3% for A 2j 3 for 5% < A 2> 3 < 25%, and therefore, within this range, 
the effect of coannihilation on the freeze-out temperature and g*(xf) can be ignored and one 
approximate the factor k by the ratio (I a (xf) + Ib(xf)v 2 )/bnv' 2 . However for A 2>3 << 5%, 
Xf can be shifted by more then a factor of two, which results in a freeze-out temperature of 
about 50 MeV. Hence, for Ni to have the observed cosmological relic density, it must must 
be that 

« E l^l 2 = (H4.04 ± 3.56 ) x ID"' x (^) 2 ^V^ , (3.9) 

ct,/3 

In Fig. we plot the ratio Vt' N h 2 /flN 1 h 2 versus m^, where Qjvi^ 2 (^jv ^ 2 ) * s the renc 
density estimated without (with) coannihilation effect. We see that for Ai < 0.05, the 
coannihilation becomes significant and leads to an increase in the relic density by more then 
50%, whereas for Ai < 0.01, this increase is almost factor of three. 

C. Indirect Detection constrains 

Before closing this section, we would like to comment on the detection of the dark matter 
in our model. 

• The direct detection: Since the interactions of N± involve only leptons, the iVi-nucleon 
scattering is absent at the tree level. Moreover, it can not scatter via electromag- 
netic interaction since the dipole moment for Majorana particles vanish identically. 
However, if Ni and N 2 (or N 3 ) are quasi-degenerate, then a transition magnetic 
dipole moment can be generated radiatively, and in that case an inelastic scatter- 
ing Ni +p — > N 2 +p is possible, provided that Ai < 10~ 6 ( 10 GcV ) . It is quite unlikely 
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FIG. 5: The ratio Q' N h 2 /Q^h 2 versus the lightest RHN mass m^, where Q^h 2 (Q' N h 2 ) is the 
N\ relic density without (with) the coannihilation effect. 

that such a tiny degeneracy will be stable under the radiative correction to m^ 1 and 
rriN 2 generated via one loop diagrams by the exchange of 5*2 and charged leptons. 

• The indirect detection: The Ni annihilation rate into leptons has a helicity suppressed 
s-wave term (i.e. oc m 2 /m 2 Ni ) and a p-wave contribution. For m^ > 100 GeV, both 
terms can be of the same order for g\ T ~ 1. Unfortunately, with the dark matter 
velocity in the galactic halo of the order 10 -3 , the cross section is too small to have a 
chance for the annihilation products of N\ to be detected. This helicity suppression 
of the s-wave will be lifted if the final state have in addition a spin one particle. 
This is the case for the internal bremsstrahlung (IB) processes, where a photon is 
emitted from the final state charged leptons (FSR) or from charged mediator in the 
t-channel propagator (VIB ) 3 . The later exhibits an enhancement for large photon 
energies if the dark matter and the particle in the propagator, whereas the FSR is 
dominated by the photons that are approximately collinear with either £ a or £p. The 
annihilation process N\Ni — > £ a £p^ not only could have a larger cross section then 
N\N\ — y £ a £/3 but also leads to a gamma ray signal with a sharp spectral features that 



is potentially observable at future gamma-ray telescopes. Very recently 



It has 



been shown that in this model, a DM with the mass mNi ~135 GeV could lead to 



3 Although the soft and collinear FSR is logarithmically enhanced, it is helicity suppressed, and thus 
typically smaller then VIB . However, they have to be both included to obtain a gauge invariant result. 
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the following three effects: (1) a wide internal bremsstrahlung bump with maximum 
of 90%mAT! ~120 GeV, (2) a 77 line around E 1 = m^ ~135 GeV, and (3) a Z7 line 
at E y = m Nl (l — m 2 z / (4m 2 Ni )) ~119.6 GeV. These features together provide a good 
fit to the gamma rays excess observed in the Fermi-LAT data 23]. 



IV. THE HIGGS DECAY CHANNELS h -> 77 AND h -> jZ 



Recently, ATLAS 24J and CMS 25[ collaborations have announced the observation of a 
scalar particle with mass ~ 125 GeV at about 5 a confidence level. The question is whether 
or not this is really the SM Higgs or some Higgs-like state with different properties. Indeed, 
the fit of the data by the ATLAS collaboration seems to show an excess in h — > 77 events by 
more than 50% with respect to the SM, while the updated CMS analysis is consistent with 
the SM. Defining i? 77 to be the decay width of h — > 77 scaled by its expected SM value, we 
find that 
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77 



1 + 



6-1 60 



(4.1) 



2 A 1 (r w ) + N c Q 2 A 1/2 (r t ) 

where Tx = ^ft/^rnJo with nix is the mass of the charged particle X running in the loop, 
N c = 3 is the color number, and Qt is the electric charge of the top quark in unit of |e|. The 
loop amplitudes Ai for spin 0, spin 1/2 and spin 1 particle contribution are given by [26| 

A 0) = -x~ 2 [x - f (x)] , 
A 1/2 (x) = 2x~ 2 [x + (x-l)f Or)] , 
A 1 (x) = -x~ 2 [2x 2 + 3x + 3 (2x - 1) / (x)] , (4.2) 
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(4.3) 



The effect on the B{h — > 77) charged scalar singlets will depend on how light are S 



1,2' 



the sign and the strength of their couplings to the SM Higgs doublet. For instance, an 
enhancement can be achieved by taking Ai and/or A2 to be negative. 

Another Higgs decay channel that could be modified due to these extra charged fields, is 
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the process h — > jZ, where similarly the effect is parameterized by 
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where (x — m2 z/^ m x-i an d the ^' s functions here are given by 26] 
A (x, y) = h (x, y) , 
4i/2 (x, y) = h {x, y) - h (x, y) , 
Ai (x, y) = [(1 + 2x) tan 2 9 W - (5 + 2x)] A (x, y) + 4 (3 - tan 2 W ) J 2 (x, y) , 

with 
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(4.6) 



(4.7) 



In Fig. [61 we present i? 77 versus i? 7 z for randomly chosen sets of parameters where the 
charged scalars are taken to be heavier than 100 GeV, the Higgs mass within the range 124 < 
rrih < 126 GeV, and the condition of a strongly first order phase transition is implemented 
(see next section). In our numerical scan, we take the model parameters relevant for the 
Higgs decay to be in the range 



A < 2, |Ai, 2 | < 3, m\ 2 < 2 TeV 2 



(4.8) 



where the Higgs mass is calculated at one-loop level. An enhancement of B(h — > 77) can be 
obtained for a large range of parameter space, whereas B(h — » jZ) is slightly reduced with 
respect to the SM. It is interesting to note that if one consider the combined ATLAS and 
CMS di-photon excess, then R 1 z is predicted to be smaller than the expected SM value by 
approximately 5%. 

V. A STRONG FIRST ORDER ELECTROWEAK PHASE TRANSITION 

It is well known that the SM has all the qualitative ingredients for electroweak baryo- 
genesis, but the amount of matter-antimatter asymmetry generated is too small. One of 
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FIG. 6: The modified Higgs decay rates B{h — > 77) vs B(h — > jZ), scaled by their SM values, due 
to the extra charged scalars, for randomly chosen sets of parameters. The magenta (yellow) line 
represents the ATLAS (CMS) recent measurements on the h—> 77 channel, while the blue one is 
their combined result. 

the reasons for that is due to the fact that the electroweak phase transition (EWPT) is 
not strongly first order, which is required to suppress the sphaleron processes in the broken 
phase. The strength of the EWPT can be improved if there are new scalar degrees of free- 
dom around the electroweak scale coupled to the SM Higgs, which is the case in the model 
that we are considering in this paper. 

The investigation of the transition dynamics and its strength requires the precise knowl- 
edge of the effective potential of the CP-even scalar fields at finite temperature |27| . The 



zero temperature one-loop Higgs effective potential is given in the DR scheme by 

v^ {h) = > - £ k . + 5>S> L (?M) _ |) , (5 ,) 



4! 2 ^ ' 64vr 2 V V A 2 

i 

where h = (y/2Re(H°) — v) is the real part of the neutral component in the doublet, rij are 
the field multiplicity, mf{h) is a field-dependent mass squared which are given in Appendix 
B, and A is the renormalization scale which we choose to be the top quark mass. At tree- 
level, the parameter /i 2 in the potential is given by fi 2 = Xv 2 , but if the one-loop corrections 
are considered, the parameter /x 2 is corrected by the counter-term 



2 ^ m dm 2 m 2 f fm 2 



v dh32n 2 \~~\A 2 J ~J h=v ^ +s / (5 ' 2) 

For instance, the one loop correction to the Higgs mass due to the charged singlets, when 
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neglecting the Higgs and gauge bosons contributions, is 

< = 2^ + Dgln=|, (5-3) 



■( 



where the first term on the right hand side of the equation is the Higgs mass at the tree 
level. If one takes ms 1 = ms 2 = 2m t and Ai = A2, then the Higgs mass is exactly 125 GeV 
for A = 10 _1 , 10" 2 , 10" 3 if Ai = 1.82, 3.68, 3.82, respectively. Note that these values are 
still within the perturbative regime. On the other hand, these extra corrections could be 
negative and may relax the large tree-level mass value of the Higgs to its experimental value 
for A large. Therefore, it is expected that these extra charged scalars will help the EWPT to 
be strongly first order by enhancing the value of the effective potential at the wrong vacuum 
at the critical temperature without suppressing the ratio v(T c )/T c , and therefore avoiding 
the severe bound on the mass of the SM Higgs. 

In order to generate a baryon asymmetry at the electroweak scale 



the anomalous 



violating B + L interactions should be switched-off inside the nucleated bubbles, which 



implies the famous condition for a strong first order phase transition 29] 



v{T c )/T c > 1, (5.4) 

where T c is the critical temperature at which the effective potential exhibits two degenerate 
minima, one at zero and the other at v(T c ). Both T c and v(T c ) are determined using the 



full effective potential at finite temperature J27I ] 



with 



and 



V eff (h,T) = V T=0 (h) + ^Y, n i J B,F {m 2 /T 2 ) + V ring (h,T); (5.5) 

i 
poo 

Jb,f (a) = / x 2 log(l =F exp(-v / x 2 + a)), (5.6) 

Jo 

V ring (h,T) = -^ J> {rhKKT) - m\{h)} , (5.7) 

i 

where the summation is performed over the scalar longitudinal gauge degrees of freedom, and 
m 2 (h, T) are their thermal masses, which are given in Appendix B. The contribution (15. 7\\ is 
obtained by the resummation of an infinite class of of infrared divergent multi-loops, known 
as the ring (or daisy) diagrams, which describes a dominant contribution of long distances 
and gives significant contribution when massless states appear in a system. It amounts 
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FIG. 7: The dependance of the Higgs vev scaled by the zero temperature value v = 246 GeV 
, on the temperature below (solid lines) and above (dashed lines) the critical temperature for two 
benchmarks, where the red (blue) one corresponds to small (large) A value and the positive (negative) 
scalar contributions in &5.3\) relax the Higgs mass to it experimental value. 

to shifting the longitudinal gauge boson and the scalar masses obtained by considering 
only the first two terms in the effective potential 30] . This shift in the thermal masses 
of longitudinal gauge bosons and not their transverse parts tends to reduce the strength 
of the phase transition. The integrals f)5.6p is often estimated in the high temperature 
approximation, however, in order to take into account the effect of all the (heavy and light) 
degrees of freedom, we evaluate them numerically. 

In the SM, the ratio v(T c )/T c is approximately (2m^ + m 3 z ) / (irvm%), and therefore the 
criterion for a strongly first phase transition is not fulfilled for m^ > 42 GeV. However, if 
the one-loop corrections in ( 15. 3ft are sizeable, then this bound could be be relaxed in such 
a way that the Higgs mass is consistent with the measured value at the LHC This might 
be possible since the extra charged singlets affect the dynamics of the SM scalar field VEV 
around the critical temperature 31]. 

This is shown in Fig. [7J where one sees the evolution of v(T) with respect to the tem- 
perature. In contrast to the SM, where the EW vev decays quickly to zero just around 
T ~ 100 GeV, here it is delayed up to TeV due to the existence of the extra charged scalars. 
This can be understood due to the fact that the value of the effective potential at the wrong 
vacuum ( < h >= 0) is temperature-dependant through the charged scalars thermal masses 
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FIG. 8: In the left figure, the critical temperature is presented versus the quantity v c /T c in j5.4\ ). 
In the right one, the relative contribution of the one-loop corrections (including the counter-terms) 
to the Higgs mass versus the parameter A. 

in the symmetric phase. The evolution of the effective potential at this (wrong) minimum 
makes the transition happening at T > 100 GeV, while the Higgs vev is slowly decaying 
with respect to the temperature as shown in FigJTl 

In Fig. [HJ we show two plots: one for v(T c )/T c versus the critical temperature, and the 
second one for the dependence of the one loop correction to the Higgs mass on its quartic 
coupling for the same sets of parameters used in Fig. [6] in the previous section. It is worth 
noting that the parameters rji, r] 2 and 7712 in (12.21) do not play a significant role in the 
dynamics of the EWPT, and therefore we fixed them in such a way to avoid the existence 
of electric charge breaking minima. 

From the left panel in Fig. [HJ we can see that one can have a strongly first order EWPT 
while the critical temperature lies around 100 GeV. The right panel shows that the one-loop 
contribution to the Higgs mass can be large compared to its tree-level value for small values 
of the self coupling A. For larger values of A, this contribution can be negative in order to 
bring the large tree-level Higgs mass down to its experimental value. Therefore, the EWPT 
can easily be strongly first order without being in conflict with the measured value of the 
Higgs mass. 

Another issue in the investigation of the EWPT that could have impacts on collider 
signatures is the possible connection between the EWPT strength and the value of the 
mass-dimension triple Higgs coupling A3 as first discussed in 32|]. In order to show this 
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FIG. 9: TTie triple Higgs coupling A3 in absolute value estimated at one loop in units of the Higgs 
vev; is shown versus the quantity v c /T c in |5.^[ j. The green line represents the tree-level value, and 
the corresponding benchmarks are the cases where different one-loop correction cancel each other. 



correlation between the EWPT strength and the enhancement on the triple Higgs coupling 
due to the non-decoupling loop effect of the additional charged singlets, we use the same 
values of the parameters of Fig. El and plot the triple Higgs coupling scaled the Higgs vev 
versus the EWPT strength, i.e., the ratio ( I5.4I) as shown in Fig. [9j 

It is clear that the triple Higgs coupling one-loop corrections could be very large with 
respect to the tree-level value for v c /T c < 2.2. 

According to the ILC physics subgroup, the triple Higgs coupling can be measured with 
about 20% accuracy or better at y/s = 500 GeV with integrated luminosity C = 500 fb^ 1 
331 ]. This implies that for large parameter space, the model can be potentially testable at 
future linear colliders. 



VI. COLLIDER PHENOMENOLOGY 

Since the RH neutrinos couple to the charged leptons, one except them to be produced 
at e~e + colliders, such as the ILC and CLIC with a collision energy \fs of few hundreds 
GeV up to TeV. If the produced pairs are the form N 2 ^N 2 ^ or NiN 2<3 , then N 2t3 will decay 
into charged lepton and S 2 , and subsequently S 2 will decay into Ni and a charged lepton. 
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If such decays occur inside the detector, then the signal will be 

fi + 2£ R , for e+e" -)• N x N 2j3 
/E + A£ R , for e+e" -> N 2 , 3 N 2>3 . 

However, for m^ > 100 GeV, it is very possible that the decay N 2 ^ — > N± + 2£ R occurs 
outside the detector, and thus escapes the detection. In this section, we assume that this is 
the case. Therefore, we analyze the production of all possible pairs of RH neutrinos, tagged 
with a photon from an initial state radiation, that is e _ e + — > NiNk'y (with i,k = 1,2,3), 
where one searches for a high pt gamma balancing the invisible RH neutrinos. 

If the emitted photon is soft or collinear, then one can uses the soft/collinear factorization 



form 34] 



- ''T/ T^ " H*,«»<D!> (e + e- -> W>) . (6.1) 

axa cos v ' 

with x = 2E 7 /y/s, here 9 is the angle between the photon and electron and a is the cross 

section ( 16. 6h evaluated at the reduced center of mass energy s = (1 — x)s. The function T 

has a universal form 

y r (x,cost^) = 5—. (6.2) 

7r x sin' 1 ^ 

Collinear photon with the incident electron or positron could be a good positive signal, 
especially if the enhancement in ( 16. ip is more significant than the SM background. 

There are two leading SM background processes: a) the neutrino counting process e _ e + — > 
vv^i from the t-channel W exchange and the s-channel Z exchange, and b) the Bhabha 
scattering with an extra photon e _ e + — > e _ e + 7, which can mimic the NiNi signature when 



the accompanying electrons or photons leave the detector through the beam pipe [35[]. In 
addition to putting the cut on the energy of the emitted photon, one can reduce further the 
mono-photon neutrino background, by polarizing the incident electron and positron beams 

such that 

K- - N„~ N e+ - N e+ 

» 50%; -^ -^ « 50%, (6.3) 



N e - + N e - ' iV e+ + iV e 

where N~ and N+ are the number densities of the left (right)-handed electrons and 

e L,R C L,R v ' 

positrons per unit time in the beam. At ^/s >> 100 GeV the process e _ e + — > vv^ is 
dominated by the ^-exchange, and hence one expect that having the electron (positron) 
beam composed mostly of polarized right handed (left handed) electrons (positron) reduces 
this background substantially, whereas the signal increases since Ni couples to the right 
handed electrons. 
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Now, let us estimate the total cross section a (e + e~ — > NiNk), which is basically the 
reverse of one of the processes which determines the effective dark matter density for coan- 
nihilation, at a collision energy of yfs. The differential cross section of e + e~ — > NiN k for the 



energy y/s is given by [181 ] 



da (e+ e ->JV a JV fc ) \9ie9%^ ft fc ( (i-Xj)(i-x k ) ^- Xi )(u-x k ) 2^x~xj \ /„ ,x 

dcostf K ik l28n s y ^y + { -_ Xb) * (t-x s )(u-x 3 ) ) ' l °^ 

with Kik = 1/2 if the two RH neutrinos are identical and equal to one if they are different, 
9 is the angle between the incoming electron and the outgoing N iy and 



x j 



m 2 Nj /s, x s = m 2 s Js, (3 ik = y (1 - x* - x k f - AxiX k 



t ii 

t = _ = 2i±a _ i (i _ p ik cos9) , u = - = x -^ - \ (1 + ft* cos 0) , (6.5) 

By integrating over cos 9, the total cross section reads 

a (e + e^ -+ iW fe ) = ^ fc ^#^^ (l + ^""TX^ 1 + !^^^In f ^ \ , 

\ IK) IK 2,2-K S I W Pik Pik W \W+l3 ik J j ' 

(6.6) 

with w = — 1 + Xi + Xfc — 2x s . In order to estimate the differential cross section of the 
process e + e~ — > NiN k ^y we integrate (16. ip over 9 taking into account the minimum value of 
electromagnetic calorimeter acceptance in the ILC to be sin^ > 0.1 36]. 

In Fig. [TU1 we show the photon spectrum for two values of collision energies y/s = 500 
GeV and 1 TeV. These plots are estimated using the factorization formula (16. ip . however, 
we obtain similar results using CalcHEP [37|. We see that for the benchmark shown in Fig. 
[TOl the heaviest RN is largely produced due to its large couplings to the electron/positron. 
Thus, for this particular benchmark the missing energy is dominated not by the DM, but 
rather by the other RH neutrinos. 

VII. CONCLUSION 

In this paper, we analyzed a radiative model for neutrino masses, generated at three loop 
level. Beside it can accommodate the neutrino oscillation data and be consistent with the 
LFV processes, it provides a DM candidate with a mass lying between few GeV up to 225 
GeV; and a relatively light charged scalar, 5^, with a mass below 245 GeV. Furthermore, we 
showed that the charged scalar singlets can give an enhancement for B (h — ?• 77), whereas 
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FIG. 10: The photon spectra from the processes e + e~ — > NiNtf) where the curves: red, green, black, 
blue, yellow, magenta correspond to (i,k)=(l,l), (1,2), (2,2), (1,3), (2,3), (3,3) respectively. Here, 
we considered the following favored mass values: m^ = 52.53 GeV, mjv 2 = 121.80 GeV, mN 3 = 
126.19 GeV, m S2 = 144.28 GeV, and the coupling values: g le = -4.19 x 10~ 2 , g 2e = 2.10 x 10" 2 
and g 3e = -6.75 x 10~ 2 . 

the decay B (h —$■ 7Z) get a small suppression, compared to the SM. We also find that 
charged scalars with masses close the electroweak scale make the electroweak phase transition 
strongly first order. Since JVi couples only to leptons, it can not be observed in experiments 
for direct dark matter searches. However it might be possible to search for such particle 
in indirect detection experiments, such as Fermi-LAT, and at future linear colliders, such 
as the international linear collider (ILC). Thus, for this particular benchmark the missing 
energy is dominated not by the DM, but rather by the other RH neutrinos. 
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Appendix A: Exact Neutrino Mass 

According to Fig. 12. 2\ the neutrino mass matrix element [a, (3] is given by 

iT ik (Q lt Q 2 , -Q 1 , -Q 2 ) ^^-Pl (~2if m ) m ^-, (Al) 

with Pl,r = (l=F75)/2, and the effective vertex Tik(pi-,P2, &i, k 2 ) is a function of the momenta 
Pi i2 and ki t 2 for charged leptons and scalars, respectively, which is given by 

ir ik (p!,p2, h, k 2 ) = J (-i 9ij ) PL % + _ m J^ Pl (-ig kj ) (fc+pi) t_ m | 2 HA.) (fc _ P2) ^_ m | 2 - 

(A2) 
Then the neutrino mass matrix element [a, 0\ is 



{M v ) a p = &fcdfpkm ii mt k \ s g ij g kj m Ni P L I {2w) i / k2 _ m2 n * 



d 4 fc 1 



d A Qi 1 / d 4 Q 2 1 fpo\ 

(2,r) 4 (Q?-mJ.)(Q?-m| l )((fc+Q 1 ) 2 -m| 2 ) J (2^) 4 (Q2_ m |J ( Q 2_ m | J ( (fc+(?2) 2_ m 2 J ' l J 



and since we have 



i 

(Q 2 -m2)(Q2- m 2) " m 



f^^ \Q 2 -mi Q'-rr^J ' ( A4 ) 



thus 

d 4 Q 1 



i (B (k 2 , m 2 , m 2 2 ) - B (fc 2 , m§, m|)j 

(^) T (0 2 - m g)(Q 2 -m?)((fc+Q) 2 -m|) 167r 2 m 2 ' l j 



where the Bn Passarino-Veltman function is 
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1 f 1 

B (k 2 ,ml,ml) = / dxln 

e Jo 



-x(l—x)k 2 +(l—x)m\+xm 2 l 
7fi ' 



(A6) 



Then, by neglecting the charged lepton masses, and making Wick rotation, we get 



(M v )ap = K^Q" f ai fpk9ij9kjF (m 2 N . /m| a , m 2 Sl /m 2 s \ , (A7) 



(47r2) J m S2 

with 



oo 



1 \ 2 



a x(l-x)r+(l-x)/3+x 



FM) = wJ dr ^\L dxln %-*»% 
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Appendix B: Thermal Masses 

The thermal masses are given as fh 2 (h) = m 2 (h) + IL; (T), where m 2 (h) and IIj (T) are 
the field-dependant masses and the thermal self-energies respectively. In this model, the 
field-dependant masses are given by 



™% = m 1 + Aii-, m 2 w s(h) = g\\, m 2 w3 _ B {h) = g 2 g x \, m\(h) = g\\, (Bl) 



where the diagonalization of the {W 3 — B} matrix gives the massless photon and Z mass 
m 2 z = (g 2 + g 2 ) h 2 /4. The thermal self-energies, that are generally estimated in the high 
temperature approximated as LTj ~ T 2 . In this model, these thermal self-energies are given 
by 

Ai A2 
~6 ~6~ 



n, 


rr ^2 / A , 3^| + g\ , y 2 

x ~ \ 12 + 16 + 4 


i w ± 


= U W 3 = —g 2 T , n^ 3 _ B = 0, 



Hb 


= —g 2 T 2 + -g 2 T 2 + -g 2 T 2 . 
6 yi 6 yi 6 yi 


n Sl 


^{M+f}. 


n S2 


T 2 f A 2 7?12 V2 1 

13 6 6 J ' 



(B2) 



The last two terms in LT/i, II X , IL3, IT^ and LTs^ represent the thermal loop contributions of 
Sf and S 1 ^ respectively. 
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